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Abstract 

The two-body scattering length a for ^He atoms is much larger than their effective range r^. As a 
consequence, low-energy few-body observables have universal characteristics that are independent 
of the interaction potential. Universality implies that, up to corrections suppressed by r^/a, all 
low-energy three-body observables are determined by a and a three-body parameter A*. We give 
simple expressions in terms of a and A^, for the trimer binding energy equation, the atom-dimer 
scattering phase shifts, and the rate for three-body recombination at threshold. We determine A* 
for several ^He potentials from the calculated binding energy of the excited state of the trimer 
and use it to obtain the universality predictions for the other low-energy observables. We also 
use the calculated values for one potential to estimate the effective range corrections for the other 
potentials. 
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I. INTRODUCTION 



The interactions of nonrelativistic particles with extremely low energies, such as cold 
atoms, are determined primarily by their S-wave scattering length a. Typically \a\ is of 
the order of the natural length scale I associated with low-energy interactions, which for 
short-range interactions is given by the range of the potential. If \a\ is much larger than I, 
however, low-energy atoms exhibit universal properties that are independent of the interac- 
tion potential. In the two-body sector, universality implies that low-energy observables are 
determined by the single parameter a up to corrections suppressed by //|a|. In particular, 
the cross section for low-energy atom-atom scattering is a simple function of E and a only. 
If a > 0, there is also a shallow two-body bound state (the dimer) whose binding energy is 
determined by a: B2 = fi^ /ma?, where m is the mass of the atoms 0]. 

Efimov showed that if the two-body scattering length is large, there are also universal 
properties in the three-body sector 2]. The most remarkable is a sequence of three-body 
bound states (trimers) with binding energies geometrically spaced in the interval between 
fi^ /ma? and f? /mt^. In addition to the binding energies of these Efimov trimers, the low- 
energy three-body observables include scattering rates for three atoms and, if a > 0, an atom 
and a dimer. The consequences of universality in the three-body sector are simplest if the 
potential supports no deep two-body bound states with binding energies of order h^/mP. 
In this case, the low-energy three-body observables are determined by a and a single three- 
body parameter up to corrections suppressed by Z/|a|. A simple physical observable that 
can be used as the three-body parameter is the binding energy of the Efimov state closest to 
threshold. Alternatively, the three-body parameter can be specified by a boundary condition 
on the three-body wave function at short distances. In Ref. js!], the authors introduced a 
more abstract three-body parameter A^, with dimensions of wave number (defined in Eq. (jiup 
below) which is particularly convenient for quantitative calculations within the effective field 
theory approach. 

A large two-body scattering length can be obtained by fine-tuning a parameter in the 
interatomic potential to bring a real or virtual two-body bound state close to the 2-atom 
threshold. The fine tuningcan be obtained experimentally by adjusting an external variable, 
such as a magnetic field |3] or an electric field Large scattering lengths for ^^Na and 
®^Rb atoms have been obtained in the laboratory by tuning an external magnetic field to the 
neighbourhood of a Feshbach resonance 0]. The fine tuning can also be provided accidentally 
by nature. A prime example is the ^He atom, whose scattering length a ~ 104 A |3] is much 
larger than its effective range ~ 7.3 A which can be taken as an estimate of the natural 
low-energy length scale /. 

The large scattering length for ^He makes this atom an ideal example of universality. 
The experimental information on low-energy ^He atoms is rather limited. Using diffraction 
of a molecular beam of small ^He clusters from a transmission grating, the bond length 
of the "^He dimer has recently been measured 0: (r) = (52 ± 4) A, which is an order 
magnitude larger than the effective range. The scattering length a = (l04t ^g]A and the 
dimer binding energy B2 = (^l.ltoj) mK were derived from the measured bond length using 

the zero range approximation 0. The "^He trimer and several larger ^He clusters have been 
observed |9lllo|. but no quantitative experimental information about their binding energies 
is available to date. However, there have been extensive theoretical calculations of the few- 
body problem for ^He using modern two-body potentials 11, 2|. Theoretical calculations of 



trimer binding energies have also improved, so that they now have several digits of accuracy 
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isi 3 3 ^ • They indicate that there are two trimers: a ground state with binding energy 
B^^ and an excited state with binding energy -83^''. The ground state binding energies of 



larger ^He clusters have been calculated using the diffusion Monte Carlo method Their 
excited state binding energies have been calculated usin g a combination of Monte Carlo 
methods and the hyperspherical adiabatic approximation [l7|. There have also been some 
calculations of three-body scattering observables. The S-wave phase shifts for atom-dimer 
scattering have been calculated in Refs. 0, 01 ■ The three-body recombination rate has been 
calculated at threshold [18. 19l Him and as a function of the energy |22l |. 

Universality implies that in the limit of large scattering length all low-energy three-body 
observables are determined by the scattering length a (or equivalently the dimer binding 
energy B2) and one three-body parameter such as A*. In order to determine this parameter, 
a low-energy three-body observable is required as input. We can take advantage of the 
accurate calculations of the trimer binding energies for modern ^He potentials by using B^^^ 
to determine A* for ^He. Once this parameter is determined, universality can be used to 
predict all other low-energy three-body observables for ^He atoms. 

We can also exploit accurate theoretical calculations for older ^He potentials to demon- 
strate the nontrivial nature of universality in the three-body sector j23l.l2^. Scaling variables 
are dimensionless combinations of physical observables. Universality implies that three-body 
scaling variables are nontrivial universal functions of aA*. By eliminating the dependence on 
A*, we can express one scaling variable as a universal function of any other scaling variable. 
The various ^He potentials have slightly different scattering lengths a and also different val- 
ues of A*. Therefore, if one scaling variable is plotted as a function of a second, the points 
for various ^He potentials should all lie along a universal line. Frederico, Tomio, Delfino, 
and Amorim calculated the scaling function relating B^^ / B^^ to B2/ Bf"* and showed that 
the points for various ''He potentials lie close to the universal scaling curve ji^ 2^ . 



In this paper, we collect all the information that is currently available on the universal 
properties in the three-body system of ^He atoms. We go beyond the work in Refs. 0, 
2^ in various respects. We give explicit parametrizations in terms of a and the three- 



body parameter A,, for many low-energy three-body observables, including Efimov binding 
energies, atom-dimer scattering phase shifts, and the three-body recombination rate. We 
also calculate additional scaling functions and estimate the effective range corrections. In 
Sectionini we briefly review universality in the two-body sector. In Section lTTTl we discuss the 
universal properties of the trimer binding energies. We determine the three-body parameter 
A* for various ''He potentials from the excited state binding energy B'i'^ and use universality 
to predict the ground state energy Bf\ We also demonstrate universality by exhibiting the 
correlation between the scaling variables Bf^ / B2 and B^^'^ / B2 for various ^He potentials. In 
Sections IIVI and IVi we discuss the universal properties in atom-dimer scattering and three- 
body recombination, respectively. We use universality to predict the S-wave scattering 
length ai2 and effective range Ts^u for atom-dimer scattering as well as the rate constant for 
three-body recombination at threshold. We also demonstrate universality by exhibiting the 
correlations between the scaling variable 012/0 and the energy scaling variables for various 
^He potentials. Section |VT1 contains a summary and the outlook for using universality as the 
basis for a quantitative description of cold atoms. 
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II. TWO-BODY SECTOR 



We begin by reviewing the universal properties in the two-body sector for ^He |pL|]. We 
use the phrase "low-energy" to refer to energies close to the threshold for free atoms. There 
is a natural length scale / for low-energy observables. For a short-range potential, / is set by 
the range. If the potential has a van der Waals tail, V{r) — > — Cg/r^, the natural low-energy 
length scale is / ~ {mC^/U^Y^^ . The natural energy scale for low-energy observables is 
h^/mf. For ^He, Cg = 1.46 a.u. which leads to Z = 5 A and Ji^/mP = 0.4 K. Universality 
occurs because a parameter in the two-body potential has been tuned such that the scattering 
length is unnaturally large. The scattering length a = 104 A f?\ for ^He is much larger than 
/. We can interpret the large scattering length as the result of an accidental fine tuning by 
nature of either a parameter in the two-body potential, such as its overall strength, or of 
the mass of the ^He atom. The ^He ^om has a mass that is only about 3/4 of that of ^He, 
and its scattering length is —7.1 A [11], close in magnitude to the natural low-energy length 
scale. 

If \a\ ^ /, universality can be used to describe the low-energy observables for atoms 
in spite of the fact that the van der Waals tail makes the potential long-range. Explicit 
expressions for the scattering length and effective range for a potential with a van der Waals 
tail have been derived in Ref . [i^ . A long-range potential introduces nonanalytic behavior in 
the dependence of the scattering amplitudes on the wavevector k. However if the potential 
falls off like 1/r®, this nonanalytic behavior enters first at 4*^ order in k. Such effects 
can not be reproduced by a short-range potential. Fortunately, their effects on low-energy 
observables are suppressed by 4 powers of I /a. We will focus on the universality predictions 
at leading order 'ml/a and also on the effective range corrections that are first order in I /a. 
At this level of accuracy, the effects of the van der Waals tail on low-energy observables 
can be reproduced by a short-range potential. Realistic interatomic potentials will therefore 
exhibit the same universal characteristics as short-range potentials. 

The two-body observables are the differential cross sections for two-body scattering and 
the binding energies for two-body bound states. The differential cross section for the elastic 
scattering of two identical spinless bosons with total energy E = Ti^k'^ /m has the general 
form 

da 



2L + 1 ^^^^^^^ 



(1) 



where ^^(fc) is the phase shift for the L'th partial wave. The total cross section is obtained 
from Eq. ((T)) by integrating over a solid angle of 27i to avoid counting identical particles 
twice. At low energies, the cross section is dominated by the L = "S-wave" term. The 
effective range expansion of /ccot5o(^) has the form 

kcot6o{k) = + ^rsP + . . . . (2) 

The first two coefficients define the scattering length a and the effective range r^. The 
natural size for these coefficients is the low-energy length scale /. The binding energies of 
the two-body bound states are determined by the poles of the scattering amplitude. For 
S-wave bound states, the binding energies are ^2"^ = ^^^n/""^; where k„ is a solution to the 
equation 

iKcot6o{iK) + K = . (3) 
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The natural energy scale for a bound state close to threshold is fi^/mP. 

If the scattering length a is unnaturally large, the differential cross section exhibits uni- 
versal behavior at energies small compared to h^/ml"^: 



da 4a 



kl < 1, lal > / . 



(4) 



The leading correction comes from the effective range. At wave number k of order (r^a)"^/^ 
or smaller, the error is of order r^/a. At larger wave numbers, the error increases like k'^r'l 
and becomes of order one at k of order l/r<.. Note that the differential cross section is 
determined as a function of k by the single parameter a. 

If a is large and positive, there is one additional low-energy observable. There is a 
shallow two-body S-wave bound state that we will refer to as "the dimer". Up to corrections 
suppressed by //a, its binding energy is 



Bo 



ma'' 



a:^l. 



(5) 



Alternatively, if we take B2 as input, universality gives a prediction for the scattering length: 

n 



(6) 



The leading correction to the universal prediction for B2 in (jSj) comes from the effective 
range. If we truncate the effective range expansion Q after the k"^ term, Eq. Q is a 
quadratic equation with two solutions: 



5?) 



2 



m rt 



(7) 



The solution B2^^ is an artifact of the truncation. We would expect a state with such an 
energy only if the higher order terms in the effective range expansion are unnaturally small. 



The solution B2 
Tg/a, we obtain 



(-) ; 



is the binding energy of the shallow dimer. If we expand to first order in 



B. 



(-) 



1 + 



(8) 



The corrections from higher orders in the effective range expansion are suppressed by 
/^/a^ and should therefore be comparable in magnitude to the r^/a^ correction. 

We now consider the two-body observables for ^He. In Table IH we give the calculated 
scattering length jl^, the effective range 0], and the dimer binding energies jl5l | for four 



commonly used potentials: two modern potentials LM2M2 and TTY jlj], and two older 
potentials HFDHE2 Q and HFD-B 

In Table HJ we also give some simple theoretical predictions for two-body observables. We 
give the prediction as for the scattering length obtained from (jH)) by using B2 as input. We 
also give two predictions for the dimer binding energy B2 using the scattering data a and 
as input. They are the universality prediction in Eq. (0) and the prediction including the 
first-order effective range correction in Eq. (jHJ. We can obtain estimates of the theoretical 
errors in approaches based on the universality at large a by comparing those approximations 
with the calculated value of B2. The universality prediction /ma^ differs from B2 by at 
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Potential 


a 




B2 


as 


ma'' 


Z 


HFDHE2 


124.65 


7.396 


0.83012 


120.83 


0.7800 


0.8263 


HFD-B 


88.50 


7.277 


1.68541 


84.80 


1.5474 


1.6746 


LM2M2 


100.23 


7.326 


1.30348 


96.43 


1.2064 


1.2946 


TTY 


100.01 


7.329 


1.30962 


96.20 


1.2117 


1.3005 



TABLE I: Two-body observables for four ^He potentials. Lengths and energies are given in A and 
mK, respectively. The first three columns show the calculated scattering lengths a 0], effective 
ranges j^, and dimer binding energies The last three columns show the universality 

predictions as for a using B2 as input, the universality predictions for B2 using a as input, and 
the predictions for B2 including the first-order effective range correction. (For ^He, the conversion 
constant is h^/m = 12.1194 KA^.) 

most 8% and the errors decrease to at most 0.7% if the first order effective range correction 
is included. This suggests that predictions of low-energy observables based on universality 
should have an accuracy of about 10% and that that one may be able to reduce the errors 
to about 1% by including effective range corrections. 



III. TRIMER BINDING ENERGIES 

The most dramatic prediction of universality in the three-body sector with large scattering 
length is the existence of Efimov states They are a sequence of shallow three-body 
bound states with binding energies much smaller than /mt^. If a parameter in the two- 
body potential is tuned such that a ±00, the number of these states increases roughly 
as ln(|a|//)/7r. The spacing of the deeper states is roughly geometric, with the ratio of 
successive binding energies approaching 515. The sug gestion that the excited state of the 



^He trimer is an Efimov state was first made in Ref. 2S|. Accurate calculations using modern 
potentials support that interpretation jl^, l^, H, 2^ • In Refs. 0,11^ S^], it was argued that 



the trimer ground state is also an Efimov state. If it is an Efimov state, universality can be 
used to predict its binding energy. We will show that the resulting predictions are within 
the errors expected from effective range corrections. We will also give a definitive criterion 
for a three-body bound state to be an Efimov state. 

Efimov derived some powerful constraints on low-energy three-body observables for sys- 
tems with large scattering length . They follow from the approximate scale- invariance at 
length scales R in the region / -C -R -C |a| together with the conservation of probability. He 
introduced polar variables H and ^ in the plane whose axes are 1/a and the energy variable 
sg\i[E)\'mE\^^'^ /fi. The angular variable ^ is 



^ = — arctan(aYmi?3//i) , a>0 



= — TT + arctan(|a|y mi?3//;,) , a<0. (9) 

Efimov showed that low-energy three-body observables are determined by a few universal 
functions of the angle E,- In particular, the binding energies of the Efimov states are solutions 
to an equation involving a single universal function A(i^) Q. Efimov's equation for the 
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binding energies reads 0, Is^] 

S3 + ^ = ^ e^-"/^" exp [A (0 /so] , (10) 
where Sq ~ 1.00624 is a transcendental number that satisfies the equation 

So cosh(7rso/2) = 8 sinh(7rso/6) . (11) 

We use a three-body parameter A* that was first introduced in Ref. through a rather 
technical definition specific to an effective field theory. Efimov's Eq. (fTUIl together with the 
explicit parametrization of A(^) given below in Eqs. ()12II14|) provides an equivalent definition 
of A*. Note that we measure from the 3-atom threshold, so > B2 for a > 0. If the 
universal function A(,^) is known, the Efimov binding energies B^ can be calculated as a 
function of a and A* by solving Eq. (jlUj) for different values of the integer n. Eq. (jiup 
has an exact discrete scaling symmetry: if there is an Efimov state with binding energy 
i?3 for the parameters a and A*, then there is also an Efimov state with binding energy 
A^i?3 for the parameters X~^a and A* if A = exp[n'7r/so] with n' an integer. Because of this 
symmetry, Eq. (fTUj) defines A* only up to multiplicative factors of exp[7r/so]. If a > 0, the 
scattering length a in Eqs. Q fTUj) can be replaced by defined in Eq. ©. The change 
in the predictions for B3 when is used instead of a can be taken as an estimate of the 
theoretical error associated with effective range corrections. 

The universal function A(^) could be determined by solving the 3-body Schrodinger 
equation for the Efimov binding energies in various potentials whose scattering lengths are 
so large that effective range corrections are negligible. It can be calculated more easily by 
using^he effective field theory of Ref. [3] in which the effective range can be set to zero. In 
Ref. y , the dependence of the binding energy on a and A* was calculated for the shallowest 
Efimov state and a > 0. In Ref. 3^1, the binding energies of the three lowest Efimov states 





■ 37r 
1- 8 


--1 

4J 


A = 


S.lOx^ - 9.682; - 


2.18, 






1- 8 


37rl 
8 -1 


: A = 


= 1.17y^ + 1.97y^ 


+ 2.12y 


-8.22, y = 7r/2 + e, 




[-^, 


SttI . 
8 -1 • 


A = 


0.25^2 + 0.282 - 


9.11, 


z = (7r + 0'exp[-l/(7r + 0' 



were calculated for both signs of a and used to extract the universal function A(,^). In 
Fig.IH we plot -{mBs/h^Aiy/'^ as a function of sgn(a)(A^,|a|) for these three branches 
of Efimov states. The binding energies for deeper Efimov states and for shallower states 
near (A^,|a|)~^/^ = can be obtained from the discrete scaling symmetry. Parametrizations 
of A(,^) in various regions for ^ were obtained by fitting the Efimov spectrum [30l |: 

(12) 
(13) 
(14) 

These parametrizations deviate from the numerical results by less than 0.013. The disconti- 
nuity at ^ — ^ = Isss than 0.016. Using Eq. (fTUj) and the parametrizations 
fll2lll4p . the full spectrum of Efimov states can be calculated as a function of a and A*. 
Eq. (jiup can also be used as an operational definition of the three-body parameter A^,, 
which was originally defined in the framework of effective field theory j^. If the binding 
energy B-^ of an Efimov state is known either from experiment or by solving the three-body 
Schrodinger equation, we can determine A,, by demanding that -B3 be a solution to Eq. (jTU)) 
for some integer n. 

A given 2-body potential is characterized by values of a and A* and corresponds to a 
vertical line in Fig.Q The dashed line shown corresponds to the LM2M2 and TTY potentials 
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-5.0 

-2.0 0.0 2.0 4.0 

sgn(a)(A.|a|)"'^^ 

FIG. 1: The energy variable — (mi?3/7i^A^)^/^ for three shallow Efimov states as a function of 
sgn(a)(A*|a|)~"^/^. The vertical dashed line corresponds to the LM2M2 and TTY potentials for 
^He. 



for ^He atoms. The intersections of this line with the binding energy curves correspond to 
the infinitely many Efimov states. The two intersections visible in the figure correspond to 
the excited state and the ground state of the ^He trimer. The third bound state predicted 
by Efimov's equation has a binding energy that is approximately 515 B^^ ^ 67 K. This is 
much larger than the natural low-energy scale /mP which is 0.4 K. This state and all the 
deeper Efimov states are therefore artifacts of the limit a ^ I. 

We can use one of the trimer binding energies as the input to determine A*. It is safer 
to use the binding energy B^^ of the excited state, because it is least affected by the high- 
energy effects that cut off the Efimov spectrum. The most accurate calculations of B^^^ have 
been obtained by solving the Faddeev equations in the hyperspherical representation jl3l |. 



in configuration space [1J|, and with hard-core boundary conditions [15[- These methods 
give results that agree to within 0.6%. The results of Ref. ^Hj for B^^ and B^^ for the 
HFDHE2, HFD-B, LM2M2, and TTY potentials are given in Table M The results for 
the binding energy B^^ of the ground state of the trimer agree well with Diffusion Monte 
Carlo calculations lal which give bP = (131.0 ± 0.7) mK for the HFD-B potential and 
(125.5 ± 0.6) mK for the TTY potential. Taking the calculated dimer binding energy B2 
as the two-body input, we determine A^, by demanding that -83^'' satisfy (jlO|) with n = 1. 
Solving Eq. (jlO|) with n = 2, we obtain the predictions for B-f^ in column 4 of Table ITTl The 
predictions are only 1-4% higher than the calculated values. If we use the calculated values 
of a as input instead of B2, we obtain the predicted values of B^^ in column 7 of Table ITTl 
These values are larger than the calculated ones by 11-21%. The difference between these 
predictions and those obtained by using B2 and B^^^ as the input gives an indication of the 
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Potential 


O 


O 




(LO) 


b!^^ (NLO) 




5^°^ (LO) 


b!^^ (NLO) 


HFDHE2 


116.7 


1.67 


1.258 


118.5 


116.7 


1.364 


129.1 


119.3 


HFD-B 


132.5 


2.74 


0.922 


137.5 


[132.5] 


1.051 


159.7 


[132.5] 


LM2M2 


125.9 


2.28 


1.033 


130.3 


126.8 


1.155 


147.4 


128.6 


TTY 


125.8 


2.28 


1.025 


129.1 


125.6 


1.147 


146.4 


127.5 



TABLE II: The trimer binding energies B^ and B^ in mK measured from the three-atom thresh- 
old for four ^He potentials. The calculated values from Ref. Jj| are in columns 1-2. The values 
of aA*, the universality predictions for ^3'^^ and the predictions for B^^ including effective range 
corrections using B2 and B^^ as input are in columns 3-5. The corresponding values using a and 
^3^^ as input are in columns 6-8. Numbers in brackets were used as input. 



size of effective range corrections. The predictions are labelled LO in Table ITT| because they 
are the universality prediction at "leading order" in r^/a. 

The numerical values of A* for the four potentials in the Table |n] are nearly the same. 
If we use ub and as the input, we obtain A* = 0.0107 A'^ for the LM2M2 and TTY 
potentials. The values for the other two potentials differ by less than 3%. If we use a and 
5f ^ as the input, we obtain A* = 0.0115 A"^ for the LM2M2 and TTY potentials. The 
values for the other two potentials differ by less than 5%. The small differences between the 
values of A* for these potentials illustrates the fact that A* tends to be insensitive to the 
parameter in the potential that is tuned to make the scattering length large. 

The availability of accurate calculations of B^^ and B^^ for various older ^He potentials 
can be used to demonstrate the nontrivial nature of universality in the three-body sector. 
Different potentials that give a large two-body scattering length should correspond to dif- 
ferent values of A*. The scaling variables B^^ / B2 and B^^ / B2 are functions of aA* only. If 
we eliminate A*, we obtain a prediction for B^^ / B2 as a universal function of B^^ /B2. A 
closely related scaling function that expresses B^^ /B^^ as a function of B2/B^^ has been 
calculated bv Frederico, Tomio, Delfino, and Amorim using the renormalized zero-range 
model jl^ U^- We have reproduced their scaling function using the solution to Efimov's 
equation In Fig.|21 our calculation of the universal scaling function relating B^^ / B2 to 
Bf^ / B2 is shown as a solid line. As A* increases, one moves along the solid line to the right. 
The data points in Fig. |21 are the results from calculations with various ^He potentials. The 
filled symbols show the results from Motovilov et al. which we used to determine A* for 
each potential, while the open symbols display the results from various other calculations 
B Q, m, Isi, 32] . The results for the LM2M2/TTY, HFDHE2, and HFD-B potentials 



are displayed by the circles, squares, and diamonds, respectively. The open triangles show 
results from Ref. jl^] for four other potentials. 

The points fall very close to the universal scaling curve, with the exception of the result 
of Ref. js^] for the LM2M2 potential, which lies well below. In that paper, the overall 
strength of the potential was varied such that B2 passed through zero. The results failed 
to exhibit the Efimov effect of an accumulation of three-body bound states at threshold 
as i?2 — > 00. The numerical accuracy of this calculation has been questioned in Ref. js^. 
All the remaining results fall along the universal scaling curve. This demonstrates that the 
dominant effect of the different potentials on the trimer binding energies can be described 
by a single parameter which we can identify with A*. 
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100.0 



200.0 



300.0 



FIG. 2: The excited state energy S3 as a function of the ground state energy B^^'. The soHd 
hne is the universal scahng curve predicted by Eq. (|1U|) . The fiUed symbols show the results from 
Motovilov et al. 11511 . while the open symbols display the results from various other calculations 
0, E3, S S, & The results for the LM2M2/TTY, HFDHE2, and HFD-B potentials are 
indicated by the circles, squares, and diamonds, respectively. The open triangles show results from 
Ref. for four other potentials. 



(0) 



It is interesting to note that calculations using only the lowest adiabatic hyperspherical 
potential [1^ or the lowest orbital angular momentum wave function give results that already 
lie near the scaling curve. Including additional adiabatic potentials or higher orbital angular 
momenta L moves the point to the right along the scaling curve until convergence is reached. 
This trend is most evident in the calculation of Ref. jl^ where the partial results for Lmax = 
0, 2 can be compared with the fully converged result with L^ax = 4 (see Table 2 and 3 in 
Ref. [13). 

Note that the most accurate points in Fig. |21 all lie systematically above the scaling 
curve by approximately the same amount. One can interpret this deviation as being due to 
effective range effects. These effects are included in potential models, but the effective range 
is set to zero in the renormalized zero range model and in the effective field theory that were 
used to calculate the scaling curve. It should be possible to account for these differences 
quantitatively by taking into account effective range corrections [s^, 3^. Since r^/a ~ 



0.07, we expect that including the effective range corrections as a first order perturbation 
would shift the scaling curve by a small amount, bringing it into better agreement with the 
calculated points. 

We can take into account the effective range corrections to B^^ approximately if we 
assume that the deviation AB^^ of the leading order universality prediction from the cal- 
culated value comes almost entirely from a correction linear in r^. The calculated result 
for one potential can then be used to estimate the effective range corrections for the oth- 
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ers. Choosing B^^ for the HFD-B potential as the input and denoting the deviations of 
B^^ from the leading order universality predictions by AB^^ , we can estimate the effective 
range corrections for any other potential by 



'ABf\ _(abP\ (r,/aB)pot 



B2 \ B2 (?"s/ctB)HFD-B 

' pot \ / HFD-B 



(15) 



Since the various ^He potentials have similar values of r^, the shift AB^^ is almost the 
same for all the potentials. The resulting predictions for B^^ are shown in column 5 of 
Table |n] The corresponding prediction using a and B^^ as the input are shown in column 
8. The predictions are labelled NLO in Table ITTl because they are approximate universality 
predictions at "next-to-leading order" in Vg/a. For each of the HFDHE2, LM2M2, and 
TTY potentials, the two NLO predictions differ by less than 3%. They also differ from the 
calculated results in column 1 by less than 3%. 

We can use the results in Table HTl to give universality predictions for B^^ for potentials 
other than HFD-B both at leading order (LO) and next-to-leading order (NLO) in the 
effective range. As our best estimate, we take the average of the two predictions obtained 
by using B2 and a as the two-body input. We take the difference to be an estimate of the 
theoretical error. The universality predictions for the TTY potential are 

LO : 5^°^ = 138 ± 17 mK (TTY) , 

NLO : B^ = 127 ± 2 mK (TTY) . (16) 

The calculated value in column 1 of Table |n] lies within the error bar for both predictions. 
Note that including effective range corrections decreases the size of the error bar by an order 
of magnitude. 

The identification of the excited state of the ^He trimer as an Efimov state is well estab- 
lished 0,[l3, 3 IE IS IE 131 • We now discuss the question of whether the ground state 



of the ^He trimer should be identified as an Efimov state. The good agreement between the 
universality prediction for B^^ and the calculated value could be fortuitous. Some authors 
have used as the criterion for an Efimov state that a sufficiently large increase in the strength 
of the two-body potential should make it unstable to decay into an atom and a dimer. In- 
creasing the strength of the two-body potential decreases the scattering length. This tends 
to move the vertical dashed line in Fig.^to the right. A sufficiently large shift in the vertical 
line will move it beyond the point where that branch of Efimov states terminates on the line 
corresponding to the dimer binding energy. According to this criterion, the excited state 
of the trimer is an Efimov state but the ground state is not. However, we argue that the 
criterion for an Efimov state should not be based on how its binding energy behaves under 
a large deformation of the strength of the two-body potential, but on how it behaves under 
arbitrary small deformations of the potential. If it is an Efimov state, any small deforma- 
tion of the two-body potential will move its binding energy along the universal scaling curve 
in Fig. El The various model potentials for '^He can be interpreted as deformations of the 
"true" ^He potential.^ The fact that the binding energies for these potentials lie along the 



The short-range part of the "true" potential can be defined by the leading order Born-Oppenheimer 
approximation. Corrections to this approximation are suppressed by rrie/m, where me is the electron 
mass. 
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universal curve is convincing evidence that the ground state of the ^He trimer is an Efimov 
state. 

If the true binding energy i?3^^ of the excited ^He trimer was measured and found to 
disagree with the calculated value using modern potentials, it would indicate that those 
potentials are not sufficiently accurate to predict low-energy three-body observables. Using 
the universality approach, there would be no need to improve the potential in order to 
predict these observables. We could simply take the measured value of -83^^ as the input 
required to determine A,,. If we did choose to improve the potential, universality implies that 
to get predictions for low-energy three-body observables with errors of order r^/a, it would 
be sufficient to introduce a two-parameter deformation of the short-distance part of the 
potential and tune both parameters simultaneously so that the potential gives the correct 
values for B2 and Alternatively, we could leave the two-body potential unchanged, 

but instead introduce an artificial short-range three-body potential and tune its strength 
in order to get the correct value for -83^^ This is essentially what is done in the effective 
field theory approach; the parameter A* is varied by adjusting the strength of a three-body 
contact interaction '3|. In the case of ^He atoms, the "true" three-body potential decreases 
the binding energy B^^ of the ground state trimer by about 0.3 mK Its effect on 

-83^'' should be much smaller, because the excited state is much larger in size. Thus the 
effects of the "true" three-body potential on low-energy three-body observables should be 
very small. However, universality implies that the dominant effect on low-energy three-body 
observables from a deformation of the two-body potential that leaves the scattering length 
fixed is equivalent to the effect of adding a three-body potential. 



IV. ATOM-DIMER SCATTERING 

The differential cross section for the elastic scattering of an atom and a dimer with wave 
numbers k in the center-of-mass frame has the form 



da 



y 2L+1 ^ 

k cot — ik 



(17) 



where Siik) is the phase shift for the L^^ partial wave. These phase shifts are real-valued 
below the dimer breakup threshold at k = (4mi?2/3/;.^)^/^. Above that threshold, they 
become complex-valued because of the inelasticity from scattering into three-atom final 
states. 

If the two-body scattering length is large, the cross section for low energies E < /ma^ 
has a universal form. For L > 1, the phase shifts are universal functions of ka only. To 
the best of our knowledge, they have not been calculated. The L = phase shift (5o(A;) is also 
universal, but it depends on aA* as well as on ka. The general structure of the dependence 
on aK^ was deduced by Efimov . For k below the breakup threshold, ka cot 5q can be 
written as 

2 

ka cot 5q = ci{ka) + C2{ka) cot[so ln(aA=,) + (j){ka)] , Q < ka < —j^ , (18) 

v3 

where ci(A;a), C2(/i;a), and (f){ka) are unknown universal functions that satisfy the constraints 
ci(2/v^) = and €2(2/ VS) = 2/^3 @. Using the effective field theory of Ref. 0, we have 
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calculated these functions. The results can be parametrized as follows: 

ci(/ca) = -0.22 + 0.39 A;V- 0.17 A;V, 
C2(fca) = 0.32 + 0.82 A;V -0.14 fcV, 

(j){ka) = 2.64: - 0.83 k^a"^ + 0.23 k^a'^ . (19) 

The atom-dimer scattering length and effective range Vg^u are defined by the low-energy 
limit of the S-wave phase shift by an equation analogous to From Eqs. (fTSl ITIJ|) . we 
obtain after the use of trigonometric identities 

ai2 = a (1.46- 2.15 tan[soln(aA,) + 0.09]) , (20) 
r,^i2 = a(l.30 - 1.64tan[soln(aA,) + 1.07] + 0.53 tan2[so ln(aA,) + 1.07]) . (21) 

The atom-dimer scattering lengths ai2 for the HFD-B, LM2M2, and TTY potentials were 
calculated in Ref. and the results are given in column 1 of Table UTTl Using the values of 
A* determined in the previous section, we can predict the atom-dimer scattering length and 
compare with the calculated values. The leading order universality predictions for ai2 and 

12 are given in columns 3, 5, 7, and 9 of Table UTTl If B2 and -83^^ are used as inputs, the 



Potential 


ai2 


OB A* 


ai2 (LO) 


ai2 (NLO) 


rs,i2 (LO) 




ai2 (LO) 


ai2 (NLO) 


rs,i2 (LO) 


HFDHE2 




1.258 


87.9 


103(5) 


278 


1.364 


65.8 


101(5) 


902 


HFD-B 


135(5) 


0.922 


120.2 


[135(5)] 


6.4 


1.051 


100.4 


[135(5)] 


18.6 


LM2M2 


131(5) 


1.033 


113.1 


128(5) 


16.0 


1.155 


92.8 


128(5) 


75 


TTY 


131(5) 


1.025 


114.5 


129(5) 


14.4 


1.147 


94.0 


129(5) 


69 



TABLE III: The atom-dimer scattering lengths 01 2 an d effective ranges rs,i2 in A for four ^He 
potentials. The calculated values of ai2 from Ref. [l5| are in column 1. The values of aA*, the 
universality predictions for ai2, the predictions for au including effective range corrections, and the 
universality predictions for rs,i2 using B2 and -63^'* as input are in columns 2-5. The corresponding 
predictions using a and as input are in columns 6-9. Numbers in brackets were used as input. 



predictions for ai2 are smaller than the calculated values by about 13%. If a and are 
used as inputs, the predictions are smaller than the calculated values by about 28%. Note 
that the predictions for 12 differ by as much as a factor of five depending on whether B2 
or a is used as the two-body input. In Fig El we show the atom-dimer scattering parameters 
ai2 and rg^u as functions of aA^:. The values of QbA^: and aA^: for the TTY potential are 
indicated by the vertical dashed and dot-dashed lines, respectively. Note that 12 is positive 
definite. It achieves a minimum value that is very close to zero near aA* = 0.94 and diverges 
at aA* = 1.64. The extracted values of aA^, for ^He are fortuitously in the interval between 
the minimum and the divergence where the effective range changes rapidly with aA*. This 
leads to a large difference in the values of 12 obtained from using a or B2 as the two-body 
input. 

In Fig. m we show the S-wave scattering phase shift So{k) for the TTY potential as a 
function of the center-of-mass energy Ecm shifted by the dimer binding energy so that the 
scattering threshold is at zero energy, 

3h^k'^ , , 

Ecm + B2 = —— , (22) 
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FIG. 3: The atom-dimer scattering length ay^ja [(a)] and effective range rs^i2/a [(b)] as a function 
of aA*. The values of as^* and aA* for the TTY potential are indicated by the vertical dashed 
and dot-dashed lines, respectively. 



where k is the wave number of the atom or the dimer. This shifted energy variable has the 
advantage that the position of the scattering threshold is independent of whether B2 or a 
is taken as the two-body input. The solid and dashed lines show the universality prediction 
with B2 and a as the two-body input, respectively. The vertical dashed and dot-dashed 
lines indicate the dimer breakup threshold for B2 and a as the two-body input, respectively. 
The filled circles show the results of Ref. which were obtained by solving the Faddeev 
equation with hard-core boundary conditions. The results are in good agreement with the 
error band defined by the solid and dashed curves. 

Universality implies that the scaling variable 012 /ag is a universal function of aA*. By 
eliminating A*, we can express universal function of 012/0^. The various "^He 

potentials should all give binding energies and scattering lengths that lie along this curve. In 
nuclear physics, a similar correlation has been observed between the spin-doublet neutron- 
deuteron scattering length ai2 and the binding energy of the triton B^. Calculations of -B3 
and ai2 using various potential models for the nucleon-nucleon interaction give results that 
cluster along a line in the ai2-By, plane called the Phillips line js^j . The observed values of S3 
and ai2 also lie on that line. Modern nucleon-nucleon potentials predict a value for i?3 that 
is about 5-10% below the measured value. Accurate values for both 012 and -B3 are obtained 
by adding a short-range three-body potential and adjusting one parameter to reproduce the 
measured triton binding energy. The Phillips line in nuclear physics arises from the large 
S-wave scattering length in both the spin-triplet (r^/a ~ 1/3) and spin-singlet (rs/|a| 1/8) 
nucleon-nucleon channels (ssl. l39|. In the case of ^He, there are two Phillips lines: one for 
the "^He trimer ground state and one for the excited state. These Phillips lines are shown in 
Figs. Efa) andlHl^b), respectively. The solid line is the universality prediction from Eqs. (jlOp 
and ()20|) . As A* increases, one moves along the solid line to the left. The data points show 
the results of Ref. Q for the LM2M2/TTY (circles) and HFD-B potentials (diamonds). 
They lie close to the Phillips lines as expected from universality. For both potentials, the 
points lie slightly above the scaling curves, consistent with small effective range corrections. 
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FIG. 4: The S-wave scattering phase shifts 5o for the TTY potential using B2 as the two-body 
input (sohd hue) and using a as the two-body input (dashed hne) as a function of the center-of- 
mass energy (with the scattering threshold defined as zero energy). The filled circles show the 
fully converged calculation of Ref. 0] with L^ax = 4. The vertical dashed and dot-dashed lines 
indicate the dimer breakup threshold for B2 and a as the two-body input, respectively. 

We can use the calculated result for ai2 for the HFD-B potential from Ref. 0| to estimate 
the effective range corrections for the other potentials. Denoting the deviation of 012 from 
the universality prediction by Aai2, we can estimate the effective range correction by 



The resulting predictions for ai2 are shown in column 4 of Table IIIII The corresponding 
predictions using a and -83^^ as the input are shown in column 8. The two NLO predictions 
agree to within 2%. For the LM2M2 and TTY potentials, they agree with the calculated 
results in column 1 to within the error bars. 

We can use the results in Table UTTl to give universality predictions for potentials other 
than HFD-B for both at leading order and next-to-leading order in the effective range. 
As the prediction and the theoretical error in ai2, we take the average and the difference 
of the predictions obtained by using B2 and a as the two-body input. The universality 
predictions for the TTY potential are 



The LO and NLO predictions for the TTY potential both agree to within errors with the 
calculated value in Table IIIII Including the effective range corrections decreases the error 
by about an order of magnitude. The error in the NLO prediction in is dominated by 




(23) 



LO : ai2 = (104±21) A 
NLO : ai2 = (129 ± 5) A 



(TTY) , 
(TTY) . 



(24) 
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FIG. 5: The Phillips line for (a) the trimer ground state and (b) excited state. The solid line 
gives the universality prediction, while the data points show the results of Ref. for the HFD-B 
(diamonds) and LM2M2/TTY potentials (circles). 



the error in the calculated value for the HFD-B potential. There is no accurate calculation 
for Ts^u for any of the ^He potentials. Thus we can only give a leading order universality 
prediction for Vg^u- Because Vs^u is positive definite and because it is so sensitive to the 
precise value of aA^,, we take the universality prediction for 12 to be the geometric mean of 
the predictions obtained by using either B2 or a as input. We take the theoretical uncertainty 
to be a multiplicative factor equal to the ratio of the two predictions. The resulting leading 
order universality prediction for the TTY potential is then 

LO: r,,i2 = (32^^21) A (TTY). (25) 

In spite of the large error bars, we can predict with confidence that rs^i2 is positive, because 
the expression (pT|) is positive definite. 



V. THREE-BODY RECOMBINATION 



Three-body recombination is the process in which two of the three incoming atoms form a 
dimer and the third atom recoils to balance energy and momentum. The rate of three-body 
recombination events per unit time and unit volume in a gas of cold atoms is proportional 
to the third power of the number density 3]: u = an^. The recombination rate constant 



a; is a complicated function of the momenta of the three incoming atoms. At threshold, all 
three momenta vanish and a reduces to a number. The total three-body recombination rate 
is the sum of the rates for all the dimers. 

If the scattering length a is large and positive, there is a shallow dimer with B2 = fi^ /ma^. 
The rate constant a for recombination into the shallow dimer at threshold must be a universal 
function of aA^. It was calculated in Ref. |2l| using the effective field theory of Ref. j^. The 
result can be parametrized as 

a = 67.1 sin^fso ln(aA,) + 0.19] . (26) 

m 
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This remarkable oscillatory dependence on ln(aA*) was previously observed in calculations 
using the hidden crossing theory and the adiabatic hyperspherical representation . In 
the hyperspherical representation, the oscillatory behavior arises from interference between 
two pathways from the incoming channel on the 2^^ adiabatic potential to the outgoing 
channel on the 1*^* adiabatic potential. Effective field theory allows the argument of the sin^ 
to be determined in terms of the same three-body parameter A* that enters atom-dimer 
scattering and the trimer binding energies. 

Using the values of A^, determined in Section IIIIl we can predict the three-body recom- 
bination constant a from Eq. ()26|) . Our predictions for a for four ^He potentials are given 



Potential 




a (LO) 




a (LO) 


HFDHE2 


1.258 


3.79 


1.363 


5.95 


HFD-B 


0.922 


0.064 


1.051 


0.37 


LM2M2 


1.033 


0.45 


1.155 


1.16 


TTY 


1.025 


0.41 


1.147 


1.11 



TABLE IV: The three-body recombination constant at threshold a in 10"^'^ cm^/s. The leading 
order predictions from universality using B2 and -63^^ (a and -63^'') as input are in column 2 (4). 
The corresponding values of aA* are given in column 1 (3). 



in Table IIVI The predictions vary by more than a factor of 2 depending on whether we 
take B2 or a as the 2-body input. This large difference arises because the value of aA^, for 
^He atoms is fortuitously close to the value near aA* = 0.83 at which the sin^ factor in (|26p 
vanishes. This is illustrated in Fig. IHl where a in units of ha'^/m is plotted as a function of 
aA*. The vertical dashed and dot-dashed lines indicate the values of asA^ and aA* for the 
TTY potential, respectively. If we use a as input instead of B2, the sin^ factor is larger by 
a factor of two. Because a is positive definite and because it is so sensitive to the precise 
value of aA*, we take the universality prediction to be the geometric mean of the predictions 
obtained by using either B2 or a as the input. We take the theoretical uncertainty to be a 
multiplicative factor equal to the ratio of the two predictions. The resulting leading order 
universality prediction for the TTY potential is 

LO : a= (O.rtli) X lO^^^ cmVs (TTY) . (27) 

There have been several previous calculations of the 3-body recombination rate at thresh- 
old for ^He. Fedichev, Reynolds, and Shlyapnikov 0| calculated the rate by making a 
semi-analytic approximation to the Faddeev equations in the hyperspherical representa- 
tion. In the limit of large scattering length, they obtained a result that depends on a only: 
a = 3.9 ha'^/m. They did not observe the oscillatory dependence of a on ln(a) predicted by 
Efimov theory, so there must have been an error in their analysis. They found that solving 
the Faddeev equations numerically for the TTY potential for ^He gave corrections of about 
10%. Inserting the value of a for the TTY potential into their analytic result, we obtain a 
prediction a = 0.6 x 10~^^ cm^/s that agrees with the universality prediction in ()27|) . This 
agreement is probably fortuitous. 

Nielsen and Macek calculated the 3-body recombination rate at threshold by applying 
hidden crossing theory to the Faddeev equations in the hyperspherical representation. They 
found that in the limit of large scattering length, a could take any value between and 
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FIG. 6: The three-body recombination constant at threshold a in units of ha^ /m as a function of 
aA^ (soUd hne). The values of a^A* and aA^ for the TTY potential are indicated by the vertical 
dashed and dot-dashed lines, respectively. 



68 ha'^/m depending on some WKB phase. This is consistent with the effective field theory 
result in ()2f)j) . For a Gaussian potential that gives the same scattering length and effective 
range as the LM2M2 potential, they obtained the prediction a = 1.1 x 10^^"^ cm^/s. They 
pointed out that the result is extremely sensitive to their WKB phase because it is close to 
the value for which a vanishes. 

Esry, Greene, and Burke ji^l calculated the 3-body recombination rate at threshold by 
solving the Schrodinger equation in hyperspherical coordinates numerically for many poten- 
tials with one or at most a few two-body bound states.^ In the case a > 0, they found that 
the recombination rates could be well approximated by an empirical formula that reduces in 
the large a limit to an expression with an oscillatory dependence on In a similar to Eq. (|26|). 
In the case of ^He, their result for the HFD-B3-FCI1 potential ^ is a = 0.12 x lO^^r 
cm^/s. A new calculation in the hyperspherical adiabatic representation that also uses the 
HFD-B3-FCI1 potential has recently been carried out i^l- This calculation includes states 



with angular momentum J > 0, so that nonzero energies can be considered. At threshold, 
it agrees with the result of Ref . [i^ . 

The large uncertainty in the universality prediction for a arises because the value of aA^ 
for the ^He potentials lies fortuitously close to the zero of Eq. (f^ . In order to improve 
on the universality prediction (j26j) within the effective field theory approach, it would be 
necessary to include effective range corrections. If there was an accurate calculation of a for 
a potential for which B^^ (or another low-energy three-body observable) is known, we could 
use that result to estimate the effective range corrections for other potentials. The only 



^ The values of = 6a given in the "Present" column of Table I in Ref. [2^ must be divided by 6 in order 
to correct for a factor of 6 error in Eq. (1) of that paper. 
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calculation of a we consider accurate enough is for the HFD-B3-FCI1 potential in Ref. 2^ . 



Unfortunately, we are not aware of any calculation of -83^^ for that potential. 



VI. SUMMARY AND OUTLOOK 

The universality approach to the three-body problem with large scattering length was 
pioneered by Efimov ^ , who emphasized the qualitative insights it provides. This approach 
implies a correlation between all the low-energy three-body observables for potentials that 
have a large scattering length. The universality approach is also useful as a quantitative tool. 
It implies that up to corrections suppressed by I /a, all low-energy three-body observables 
are determined by a and a single three-body parameter. A convenient choice for the three- 
body parameter is the parameter A^, introduced in Ref. , because the dependence of some 
observables on A^. is known analytically. 

In order to determine A,,,, one three-body observable is required as input. A convenient 
choice is the binding energy of the shallowest Efimov state. Once A* is determined, all other 
low-energy three-body observables can be predicted. We used calculations of the binding 
energy B^^ of the excited state of the trimer to determine A* for various ^He potentials. 
We then used universality to calculate the binding energy Bf^ of the ground state of the 
trimer, the atom-dimer phase shifts below the dimer breakup threshold, and the three-body 
recombination constant at threshold a. We gave explicit expressions for the three-body 
recombination constant a in Eq. ()26|) and for the S-wave atom-dimer phase shifts below the 
breakup threshold in Eq. ()18p. We also gave an explicit parametrization for the universal 
function A(^) that appears in Efimov's equation ()1U|) for the trimer binding energies. 

The leading corrections to the universality predictions come from effective range correc- 
tions. If these corrections are included, there should be a systematic improvement in the 
accuracy of the predictions for all low-energy observables with errors decreasing to second 
order in I /a. The effective range corrections have not yet been calculated for the case of 
three identical spinless bosons with large scattering length. We therefore used accurate 
calculations of Bf"* and ai2 for the HFD-B potential as input to deduce the approximate 
effective range corrections in these observables for the other potentials. The resulting the- 
oretical errors are smaller than those for the leading-order universality prediction by about 
an order of magnitude as expected. Comparing with the calculated values of Bf^ and ai2, 
we see that the theoretical errors correctly reflect the accuracy of the LO and NLO univer- 
sality predictions. An actual calculation of the effective range corrections for the three-body 
observables would eliminate the need for using calculations of three-body observables for 
one potential as additional inputs. 

The leading-order universality predictions presented in this paper were obtained using 
the effective field theory of Ref. |3| , which is a particularly convenient implementation of the 
universality approach for three-body systems. More generally, effective field theory provides 
a framework for the model independent description of low-energy phenomena by exploiting 
a separation of scales in the system Using renormalization, all short- distance effects 
are systematically absorbed into a few low-energy parameters such as the scattering length 
a and A,,. As a consequence, the dependence on the relevant low-energy parameters is 
explicit, while irrelevant details of how their numerical values arise from complicated short- 
distance dynamics are omitted. Effective field theory allows for systematically improvable 
calculations of low-energy observables with well-defined error estimates. This method has 
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many applications ranging from particle physics over nuclear physics to condensed matter 
physics j4i.l43ll43. 



An important open question is how universality is manifested in the four-body prob- 
lem. Low-energy four-body observables must depend on the two-body parameter a and the 
three-body parameter A*. Are any new four-body parameters required to calculate low- 
energy four-body observables up to corrections suppressed by l/\a\l There are theoretical 
arguments in support of both answers to this question. There is a renormalization argument 
for zero-range 5-function two-body potentials that indicates that a new four-body parameter 
is required to calculate four-body binding energies |45| . On the other hand, a power counting 
argument within the effective field theory framework suggests that a four-body parameter 
should not be necessary to calculate four-body observables to leading order in //|a| [i^. This 
raises the exciting possibility of calculating the binding energy of the ^He tetramers close 
to the four-atom threshold to about 10% accuracy using a and b'^'' (or another low-energy 
three-body observable) as the only inputs. There is some circumstantial evidence in favor of 
this possibility from the four-body problem in nuclear physics. There is a correlation called 
the "Tjon line" between the binding energy of the triton (^H nucleus) and the binding 
energy B^ of the alpha particle (^He nucleus) ji^. Calculations of these binding energies 
using modern nucleon-nucleon interaction potentials give results that underestimate both 
binding energies but cluster along a line in the -B3--B4 plane. By adding a three-body po- 
tential whose strength is adjusted to get the correct value for B^, one also gets an accurate 
result for i?4 (see Ref. ji^ for some recent calculations with modern nuclear forces). 

The results for '^He presented in the paper apply equally well to other bosonic atoms 
with large scattering length as long as the effects of deep two-body bound states on low- 
energy observables are sufficiently small. By definition, a deep bound state has a binding 
energy of order 'h^/ml'^ or larger. If a < 0, any dimer is deeply bound. If a > 0, any 
dimer other than the shallow dimer with B2 ~ Ti' jmc? is deeply bound. One qualitative 
effect of the deep two-body states is that the Efimov states become resonances, because 
they can decay into an atom and a deeply bound dimer. Deep two-body bound states also 
provide additional channels for three-body recombination. Their effects can be particularly 
dramatic for a < if there is an Efimov state near the three-atom threshold, because it gives 
a resonant enhancement of the three-body recombination rate into deep two-body bound 
states |2ol . 4^. The existence of deep two-body bound states does not affect the universality 



prediction for low-energy observables in the two-body sector. However in the three-body 
sector, it implies that a third parameter in addition to a and A* is required to predict low- 
energy observables up to corrections suppressed by //|a| j^. This parameter takes into 
account the cumulative effects of all the deep two-body bound states. The modification of 
Efimov's equation for the binding energies was given in Ref. j30l |. 

The universality approach discussed in this paper is not limited to identical bosons. It 
can be applied to any three-particle system for which at least two of the three pairs have 
a large scattering length. Some examples are given in a recent review article [sOl- The 
universality predictions will depend on the three pair-wise scattering lengths, the ratios of 
the masses, and the three-body parameter A*. 

An especially promising application of the universality approach is to cold atoms in the 
vicinity of a Feshbach resonance, where the effective scattering length can be controlled by 
an external magnetic field 0]. This creates the exciting possibility of testing the unique 
oscillatory dependence of low-energy three-body observables on the scattering length that 
is predicted by universality. In this paper, we have exploited the fact that the various 
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^He potentials span a small interval of aA^,. Using a Feshbach resonance to control the 
scattering length, it might be possible to scan through an entire cycle of the oscillatory 
behavior. Among the dramatic effects that one may be able to observe are the divergence 
in ai2 near aA* = 4.3 and the zero of a near aA^, = 0.83. 

The behavior of cold atoms near a Feshbach resonance is in general a complicated coupled- 
channel problem involving the various hyperfine states of the atoms. However, from the point 
of view of universality, the coupled-channel effects introduce no additional complications. If 
one is sufficiently close to the resonance and if the energy relative to the threshold for one 
hyperfine state is small compared to the hyperfine splittings, only that hyperfine state needs 
to be included explicitly. The coupled-channel effects can be taken into account through 
the values of the low-energy parameters a, which diverges at the resonance, and A*, which 
varies slowly in the neighborhood of the resonance. 

The behavior of a Bose-Einstein condensate of atoms with large scattering lengths has 
been studied experimentally by using Feshbach resonances to tune the scattering lengths of 
alkali atoms In the low-density limit na^ -C 1, the nontrivial aspects of universality in 
the three-body sector are reflected in a small oscillatory dependence of the energy density 
of the condensate on ln(A*n^/^) j5l|. There is a possibihty that these three-body effects 
would allow the existence of stable homogeneous condensates with large negative scattering 
length [H^l- The intermediate density region where na^ ~ 1 but nl^ -C 1 is a more difficult 
problem. It is an open question whether a condensate in this region has universal properties 
that are determined by constants such as a and A^, that describe the low-energy properties 
in the few-body sectors. If there are, it may be possible to use universality to predict in 
detail the behavior of a Bose-Einstein condensate of atoms near a Feshbach resonance. 
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